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Abstract. We consider the impulsive limit of a class of space-times generalizing pp- 
waves that are given in the form M — N x where {N, h) is a Riemannian manifold 

■ of arbitrary dimension and M carries the line element ds^ = dh'^ + 2dudv + f {x)5{u)dv? 

with dh? the line element of N and 5 the Dirac measure. We prove a completeness 

■ result for these space-times. 



1. Introduction 

>: 

^ ■ Plane-fronted gravitational waves with parallel rays — pp-waves, for short — are defined 

\0 • by the existence of a covariantly constant null vector field k and are usually associated 

O, ds^ = 2dudv + dx\ + dx\ + H{x^ ^x^ ^u)du^ (1) 

on M^. These space-times model gravitational or electromagnetic waves and other forms 
of null matter and have been extensively studied (see e.g. |GP09i Ch. 17] and the 
^ \ literature cited therein). The geodesic null congruence with tangent k is non-expanding, 

\ shear-free, and twist-free and the latter property implies the existence of a family of 2- 

surfaces perpendicular to k which are interpreted as wave surfaces. Moreover, since k^.^ 
vanishes, they are planar and rays orthogonal to them are parallel. 

It should be noted, however, that Brinkmann who studied these geometries in the 
context of conformal mappings of Einstein spaces ( |Bri25] ). also included a rotational 
term (rediscovered by Bonner |Bon70] and recently studied further under the name 
gyraton |Fro07j ). as well as allowing for a general wave surface. Including the latter 
effect, i.e., allowing for a Riemannian manifold of arbitrary dimension as the wave surface 
we arrive at the following geometry {M,g): Let {N,h) be a connected Riemannian 
manifold of dimension n, set M = x and equip M with the line element 

ds"^ = dh^ + 2dudv + H{x, u)du^ . (2) 

Here dh? denotes the line element of {N, h) and u, v are global null-coordinates on the 
2-dimensional Minkowski space R^. Finally ifiMxM— T-Misa smooth function. 



On the completeness of impulsive gravitational waves 



2 



These models have been studied in a series of papers by J. Flores and M. Sanchez 
in part together with A. Candela f jCFS03l [FS031 ICFS041 [FS06] ) mainly focusing 
on causality and geodesies. These geometries allow to shed some light on some 
of the peculiar causal properties especially of plane waves (i.e., pp- waves ([1]) with 
H{x^,x'^,u) = hij{u)x^x^), see e.g. |BEE96l Ch. 13]. They turn out to be caused by the 
high degree of symmetries of plane waves and the fact that the wave surfaces of ([1]) are 
fiat M2. 

In |CFS03] space-times of the form ([2]) have been called (general) plane-fronted 
waves (PFW). However, by the geometric interpretation given above and by the analogy 
with pp-waves it seems more natural to us to call the space-times ([2]) N -fronted waves 
with parallel rays (NPW), which we shall do from now on. 

It turns out that the behaviour of H at spatial infinity, i.e., for ,, large x" is decisive 
for many of the global properties of NPWs. In order to formulate precise statements we 
recall that one says that H behaves subquadratically at spatial infinity if there exist a 
fixed point x E N, continuous functions < Ri, R2 and a continuous function p < 2 
such that 

H{x, u) < Ri{u)dP^''\x, x) + R2{u) y{x, u)eNxR. (3) 

Here d denotes the Riemannian distance function on N. Similarly we say that H behaves 
at most quadratically resp. superquadratically if p < 2 resp. p > 2. In |FS03j it has 
been shown that the causality of NPWs depends crucially on the exponent ]? in (jS]), 
with p = 2 being the critical case. In particular, NPWs are causal but not necessarily 
distinguishing, they are strongly causal if —H behaves at most quadratically at spatial 
infinity and they are globally hyperbolic if —H is subquadratic and is complete. 
Similarly the global behaviour of geodesies in NPWs is governed by the behaviour of H 
at spatial infinity. From the explicit form of the geodesic equations it follows ( |CFS03 
Thm. 3.2]) that a NPW is complete iff is complete and 

Df''x = ^V,H{x,s) 

has complete trajectories. Here D^^^ is the induced covariant derivative on and V^, 
denotes the spatial gradient. Applying classical results on complete vector fields (e.g. 
|AMR88l Thm. 3.7.15]) completeness of M follows for autonomous H (i.e., independent 
of u) in case H grows at most quadratic at spatial infinity. Clearly this implies 
completeness for at most quadratic sandwich waves, i.e., H compactly supported in 
u. 

In this work we consider impulsive NPWs (INPWs), i.e., we set H{x, u) = f{x)5{u) 
in ([2]), where 6{u) is the Dirac measure on the hypersurface {u = 0}. Impulsive pp- 
waves (for a summary see [GPOOj Ch. 20]) have been introduced by Penrose using a 
,, scissors- and-paste method" (e.g. jPen72j ) gluing two halves of Minkowski space along 
the null hypersurface {u = 0} with a warp. On the other hand, impulsive pp-waves 
arise as ultrarelativistic limits of Kerr-Newman black holes, the prototype being the 
Aichelburg-Sexl geometry ( |AS71] ). 
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The distributional term in the metric of impulsive pp-waves and INPWs makes it 
a delicate matter to mathematically deal with these space-times; for a general account 
on distributional geometries in GR see |SV06] . Therefore impulsive pp-waves have been 
treated using the nonhnear distributional geometry ( |GKOS0"T] ) built upon algebras of 
generalized functions ( |Col85] ). In particular, geodesies in impulsive pp-wave space- 
times have been considered in |Bal97t ISte98] . and in |KS99] . where an existence and 
uniqueness result for the geodesic equations has been proved. From a global point of 
view these results imply that impulsive pp-waves are geodesically complete. 

In this short note we prove a completeness result for INPWs with complete A^. 
We do so without using any theory of nonlinear distributions leaving a detailed study 
of INPWs as distributional geometries to a subsequent paper. More precisely, we view 
INPWs as geometries with a small but finitely extended impulse: Let 5^ be some smooth 
approximation of the Dirac-delta (i.e., 5^ ^ 5 weakly as e — )■ 0) and for fixed e > 
consider the metric 

dsl = dh'^ + 2dudv + f{x)5e{u)du'^ (4) 

on M, where / is an arbitrary smooth function on A^. We will show that for any geodesic 
7 in M there is eo small enough, such that 7 can be defined for all values of an affine 
parameter provided e < eo. Moreover the size of eo for which the geodesic becomes 
complete can be explicitly estimated in terms of (derivatives of) / and the initial data 
of 7. Finally, we also show that the globally defined geodesies converge to the geodesies 
of the background x which, however, have to be joined with a suitable warp at 
the shock hypersurface. 



2. The geodesic equations for INPWs 

In this section we derive the geodesic equations for INPWs and fix some notation 
to be used in the remainder of this work. We start by making precise the class of 
regularizations we use for the Dirac delta. We set / := (0, 1]. 

Definition 2.1. A net ((Je)ee/ of smooth functions on M is called a strict delta net if it 
satisfies the following three properties. 

(i) The supports shrink to zero, supp(5e) — > {0} for e \ 0. 

(ii) The integrals converge to 1, / 5e{x)(ix — )■ 1 for e \ 0. 

Jr 

(Hi) The -norms are uniformly hounded, 3K > : / |5e(2;)|da; < K Ve G /. 

Jr 

Observe that this is a very general class of approximations of 6. (Even although 
smoothness excludes "boxes", nets arbitrary close to "boxes" and even discontinuous 
regularizations are practically included by the fact that is dense in L^.) Without 
loss of generality we will always assume that supp((5e) C (— e, e) for all e G /. 
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Now let M = X be an INPW with a connected n-dimensional and complete 
Riemannian manifold and let M be endowed with the line element (jl]), where {5^e is a 
strict delta net. 

Denoting the Christoffel symbols of the Riemannian manifold (A^, h) by V^^^ 
one obtains the non-vanishing Christoffel symbols for M w.r.t. a coordinate system 
(x^, . . . , x") of N and (m, v) null-coordinates of 



^% =^^''\ foralll<z,j,A;<n, 



Kj=^)u = TTTTT'^e foralll<j <n, 



idi 

2 dx^ 

=- f6 



2 

pfc _ r^km J r 



2 ax"" 

Since all Christoffel symbols of the form T^j^ vanish we may use u as an affine 
parameter for the geodesies which hence obey the following set of n + 1 equations 

4f A. = \ V./(x.) 5,. (6) 

Here D^^^ and V^; denote the covariant derivative resp. the gradient w.r.t. h. First 
observe that equation ([5]) can be integrated once the second equation has been solved. So 
we have to concentrate on equation ([6]), which is just the perturbed geodesic equation on 
with potential / and the non-autonomous term 5^- Moreover, since the latter vanishes 
for \u\ > e the x-component of the geodesies on M will for large u coincide with the 
(unperturbed) geodesies on A^. By completeness of A^ the question of completeness 
of M reduces to the question whether all perturbed geodesies on A^ that enter the 
regularization strip a.t u = —e also leave it at u = e, that is whether the perturbed 
geodesies blow up before m = e or not. 

Bearing this in mind we apply the following procedure to solve the geodesic equation 
on M as well as to address the problem of geodesic completeness of M. We fix e > 
and impose initial data xq G A^, xq ^ T^^N "long before" the shock at n = — 1 and then 
follow the unperturbed Riemannian geodesic on A^ with this data, i.e., the solution of 

D^^^x = 0, x(— 1) = xo, i(— 1) = xo, 

which we denote by x[xo, Xq ]. By completeness of A^ this geodesic x[xo , Xq 

1 will reach 

the shock region at u = —e and until then it will also be a solution of the perturbed 
geodesic equation ([6]) with the same data, which we will denote by Xe[xo, xo]- With this 
notation we have Xe[xo,Xo] = ^^[a^Oji^o] on ] — oo, — e] and to continue xJxo,xo] into the 
shock region \u\ < e we consider the initial value problem 

(E]) with data Xe(— e) = x[xo, xo](— e), e) = x[xo, xo](— e)- 
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To prove that Xe[xo,io] extends to all values of the parameter u we only have to show 
that the latter initial value problem possesses a solution denoted by xJ\Xq, Xq] until -u = e, 
since for u> e the right hand side of (EI) vanishes and we are solving the (unperturbed) 
geodesic equation in the complete manifold A^. That is, we only have to show that no 
blow-up occurs within the shock region \u\ < e, which, in fact, will be done in the next 
section (at least for e small enough). In total we will then have the global perturbed 
geodesic x^[xo, Xq] 

r ■ 1/ \ J ^FOj^^o u < —e . . 

Xe[xo,xo\[u) = I (7) 
1^ x\xq,xq\ -e < u. 

Finally, as observed above, once we have such a solution of the x-component of 
the geodesic in M the equation for v can be integrated to give a solution for all m e M 
and we will use the following notation: for initial conditions Vq^Vq G M we denote by 
f [vo,i'o] the straight line v\vi^^v^{u) = f o + 't'o(l + i-e., a solution of ([5]) for m < — e 
and similarly ^Jvc^o] denotes a solution of (|5]) with v^i^^v^{—t) = v[vo,VQ]{—e) and 

3. Completeness 

We now show that for any geodesic in M we can choose e sufficiently small such that 
the geodesic can be extended through the shock. More precisely we prove that (using 
the notation introduced above) the initial value problem 

-^if^^e = ^'^xf{xe)Se, x^{-e) = x[xo,xo]{-e), x^{-e) = x[xo,Xo]{-e) (8) 

has a local solution defined up to it = e provided e is small enough. 

Proposition 3.1. For all xq & N, xq & T^^N there exists eo such that the initial value 
problem ^ has a solution Xf[xQ, xq\ defined up to u = e provided e < eo- 

The proof heavily rests on a fixed point argument which we provide in detail in 
Lemma [A. 2 1 in the Appendix. Here we only observe that this argument indeed provides 
the assertion of the Proposition. 

Proof: We invoke Lemma with 6 > 0, c > 0, Fi{y,zY ■= —T'^j^\y) z"- (to 
express D'^^^ in local coordinates) and F2{yY '■= \h^™'{y)^^{y) which is just \V xf in 
coordinates. Clearly Fi and F2 are smooth since / and h (and hence the Christoffel 
symbols) are assumed to be smooth. Hence Lemma IA.2I guarantees existence of a 
solution Xe[xo,Xo] of dH]) until u = a — e. So choosing eo = f , the solution 5;e[xo,xo] 
exists at least until u = +e, provided e < eg. □ 

Observe that Lemma [A. 21 also implies that the solution Xe[a^O)3:o] together with its 
first derivative a;e[a;o,xo] is uniformly bounded (in e on [— e, e]). Moreover flA.2p gives an 
upper bound on eo in terms of the initial velocity Xq and of the Christoffel symbols on 
as well as of Vxf on a neighborhood of the data xq. Next we state our completeness 
result. 
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Theorem 3.2. For all xq G N, xq G T^^N and all Vq, -Uq G M there exists eo such 
that the solution {x^[xq,Xo],v^[vo,Vq]) of the geodesic equation (01^ with initial data 
— 1) = xq, x^{—1) = xq, 1) = vo, Ve{—1) = vq is dcfincd for all G M, provided 
e < eo- 

Proof: Given zq, xq, Proposition 13.11 provides us with eo such that the solution 
of dH]) is defined for u G [— e, e] for all e < eo- In this case we hence may define 
Xe[a;o,io] as in ([7]) for all m G M so it only remains to integrate ([5]) twice to obtain a 
globally defined solution v^- Hence in total we obtain a unique globally defined geodesic 
M 9 M {x,[xq,Xo\{u),v,[vo,Vq\{u)). □ 

We point out that a in the proof of Proposition 13.11 and hence eo for which the 
geodesic is defined on all of R depends on the choice of the initial data Xq and Xq. 
Hence we can not, in general, obtain a global bound eo such that for fixed e < eo the 
manifold M is geodesically complete. There are, however, two special cases where we 
actually obtain geodesic completeness of M for e sufficiently small. First assume that 
is compact. Then we obtain a globally defined eo since xq varies in a compact set 
only and upon reparametrization we may achieve that |xo| = 1. On the other hand, if 
— / behaves subquadratically (cf. ([3])) then by the compactness of the support of 4 we 
may apply the results of |FS03] mentioned in the introduction to obtain completeness 
without even the need to invoke the fixed point argument. 

However, one may say that "in the limit e — )■ 0" we obtain a geodesically complete 
manifold, hence one may say that IPNWs are geodesically complete irrespectively of the 
behaviour of the profile function /. This is in sharp contrast to the case of extended 
NPWs where completeness depends crucially on the behavior of H at "spatial infinity" : 
the role of the x-asymptotics of H becomes irrelevant in the impulsive limit. 

However, the precise meaning of the completeness statement (i.e. the dependence 
of eo on the data) is encoded in the formulation of our theorem above. A more straight 
forward completeness result for IPNWs can be provided using nonlinear distributional 



geometry ( |GKOS0"T| IKS02] ) in the sense of J.F. Colombeau ( |Col85] ). and we will 
address this topic in a subsequent paper. 



In this section we compute the limits of the global geodesies derived above as e — > 0. We 
start by analyzing the x-component and introduce some more notation in the same spirit 
as at the end of section |2j To denote the prospective limit of xJ[Xq^ Xq\ we have to paste 
together the solution x[xo, xo] of the unperturbed equation for m < with an appropriate 
solution of the unperturbed equation for m > 0. To this end denote by x[a;o,Xo] the 
solution of the (unperturbed) geodesic equation on with data £(0) = x[xq, io](0) and 
x(0) = i;[xo,io](0) + |Vi./(a;[xo, io](0)). Finally denote the prospective limit by 



4. Limits 




(9) 
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Observe that y[xo,Xo] is a continuous curve M — > which is piece-wise smooth with a 
single break point at = 0. Moreover it is not differentiable (in general) as we have 

\m\y[xQ,XQ]{u) = x[xo,xo]{0) + ]-V^f{x[xo,xo]iO)). 

u\0 Z 



For simplicity we write Fi{y,zY '■= —T'l^^\y)z^z^ and F2 = |V^/ = ^h^''-Pi 



ij vy;^ ^ X 2 2 ^ 2'" dx' 



as 



in the proof of Theorem 13.11 and start with an auxiliary result needed throughout the 
remainder of this section. 

Lemma 4.1. The global solution x^[xq,xq] of ^ (defined in ^) satisfies 

Xe[xo, XoKeu) — )■ y[xo,XQ]{0) = x[xq,Xo]{0) uniformly on [—1,1] as e \ 0. 

Proof: To keep the notation transparent we abbreviate a;[a:o,Xo] by x and xJxcXq] 
by Xe- We have 

sup \x^{eu) — x{0)\ < sup \x^{eu) — x{eu)\ + sup |a;(eM) — a;(0)| . 

ue[-i,i] Me[-i,i]^ V ■' ue[-i,i] 

=■■(*) 

The second term goes to zero as e \ since x is uniformly continuous on compact sets. 
To estimate the first term we integrate the differential equations for x^ and x (see also 
the Appendix) to obtain 

/ \Fi{x^{r),x^{r)) — Fi{x{r),x{r))\drds + / / \F2{x^{r))\\5^{r)\drds 

< Ce^ + C\\6,\\Lie < Ce (e \ 0) , 
where we have used that by Lemma [A. 21 and are bounded independently of e. □ 

Proposition 4.2. The global solution x^[xo,Xq] of ^ (defined in Ij^) satisfies 

Xe[xo,xo] — )■ y[xo,XQ] uniformly on compact subsets o/M, 
Xe[a;o,2;o] — !■ y[xo,Xo] uniformly on compact subsets o/]R\{0}. 

Proof: Again we write x for x[xo, xq] and x^ for Xe[xo, xq] and similarly y for y[xo, xq] 
and for 5;f[xo,io]- Without loss of generality we only consider the compact interval 
[—1, 1]. We distinguish three cases: —l<u< — e, —e<u<e and e < u < 1. 

In the first case x^ = x = y on [—1, — e] (and hence x^ = x on the same interval), 
since x^ and x solve the same initial value problem. If — e < m < e the result for 
follows immediately from Lemma S]T] while for the derivative x^ there is nothing to prove 
in this case. 

Finally, for e < -u < 1 we observe that x^ = x^ and y = x solve the same differential 
equation but now with different initial conditions, namely (e), and x{e) and 

x{e), respectively. By continuous dependence on the initial data we obtain 

m8ix{\x^{u) — x{u)\^ |5;e('w) — 2;(m)|) < max(|a;e(e) — x(e)|, \xj^e) — x(e)|) c^, 
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where L is a Lipschitz constant of Fi on the compact image of [0, 1] under x, 5^, 
and it suffices to estimate the difference of the data. Indeed we have 

|xe(e)-5(e)| < |5,(e)-£(0)| + |x(0)-5(e)| 0, 

since the first term converges to zero by Lemma 14.11 and the second by continuity. 
Similarly we have 

|ie(e)-£(e)| < |£,(e)-i(0)| + |£(0)-i(e)|, 

where again the second term on the r.h.s. goes to zero by continuity. To estimate the 
first term we plug in the integral representation of x^ to obtain 

|ie(e) - £(0)1 < |5,(e) - i(0) - F2(x(0))| 

<\'x,{-e) - x{Q)\+ I \Fi{5:,{s),'x,{s)\ds 



+ 



F2{ie{s)) 5,{s) ds-F2{x{Q)) 



Now the first term on the r.h.s. vanishes in the limit since x^{—e) = x{—e) — )■ x{0). The 
second term goes to zero again by the uniform boundedness of x^ and x^- To obtain the 
same conclusion for the third term we again take into account the uniform boundedness 
of x^ and the fact that {5e)e is a strict delta net. □ 

Next we turn to the t;-component and recall that {u,v) G M.f and so we may work 
distributionally. 

Proposition 4.3. The global solution v^[vq,vq] of ^ satisfies 

^Jt;o,{;o] ^ ^;K,^'o] - ^ /(x(0))i7 - (i^(0) + ^V./^(x(0)))D,/(x(0))ii+, 

where u^{u) = uH{u) denotes the kink function and we again have abbreviated X^Xq^ Xq] 
by X. 

Proof: In addition to the abbreviations x and x^ used already above we write v for 
vft'c'i'o] and Ve for Ve[vo,vo]- Since we have Ve{u) = vq + 7)0(1 + u) + H * H * Ve{u) 
and since convolution is a separately continuous operation, it suffices to calculate the 
distributional limit of i;^. Inserting the integral representation of xi into equation ([5]) 
we obtain 



Ve{u) = - Djf{x^{u)) 6eiu) x^(-e) -- Djf{x^{u)) 6eiu) I F({xeis),x^{s)) ds 



(/) V 

in) 

- 1 D,f{x,{u)) 5,{u) Fi{x,{s))6,{s) ds ^{f{x,{u))5,{u)) . 

s ^ / ^ V ^ 
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It is easily seen that (/) ^ {Q)Djf{x{'d))5 and that (/V") f{x{Q))5 in P'(M). On 
the other hand {II) — ?■ in D'(]R) since for all test functions cf) G I?(M) we have (again 
using the uniform boundedness of x^ and x^) 



(p{u) Djf{x^{u)) 5e{u) / F({xe{r), Xe{r))drdu 



< 2eC 



Finally, we show that (///) — )■ ^Djf{x{0))F2{x{0))S. Indeed, we have 



(Piu)D,f{x,iu)) 6,iu) / Fiix,{r))6,{r)drdu - -0(O)D,/(x(O))F2^(x(O)) 



< 

+ 

+ 
< C 

+c 

+c 



(j){u)5,{u) / Fi{x,{r))6,{r)dr[D,f{x,{u)) - D,f{x{0)))du 



(t){u)5^{u) 



Fii^eir)) - Fi{x{0)))6,{r)drdu 



6e{r)drdu — -^(O) 



l^./(^(0))| 



|D,/(x(0))||F|(a;(0))| 



oo sup 

ue[-i,i] 



oo sup 

«e[-i,i] 



D,f{x,{eu) - D,f{xm 
Fi{x,{eu)) - Fi{x{0)) 



(j){u)6,{u) I 5,{r)drdu - -(f){0) 



Now the first and the second term converge to zero, again by Lemma 14.11 Finally, the 
integral term in the last line converges to zero by an elementary calculation. □ 

Summing up we have shown that the x-component of the limit is continuous but 
has a kink at the shock hypersurface. The ti-component, however, is not even continuous 
but has a jump at the shock in addition to a kink. The parameters of the kinks and the 
jump are given in terms of the profile function / and its derivatives at the point where 
the geodesic hits the shock hypersurface. So globally the geodesies on M are given by 
geodesies on the background N x M.^, which have to be joined suitably at the shock 
hypersurface. 
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Appendix 

In this appendix we detail the fixed point argument used in the proof of our main 
result. Our argument is built on a slightly sharper version of the Banach fixed point 
theorem (see |Wei52] ). Indeed the integral operator A^, used below to solve the initial 
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value problem is not a contraction on the naturally chosen Banach space and so the 
Banach fixed point theorem does not apply. 

Theorem A.l. (Weissinger's fixed point theorem) Let X be a nonempty closed subset 
of a Banach space {E, ||.||). Moreover let Yl'^=i be a convergent series of positive real 
numbers {an)n o,nd A : X ^ X a map with the property that 

\\A''{u)-A''{v)\\<an\\u-v\\ yu,veXyneN. (A.l) 

Then A has a unique fixed point. 

We may now state and proof our main technical result. 

Lemma A.2. Let Fi G C°°(M2",M"), F2 G C:~(M",M"), let xo,xo G W, let b > 0, 
c > be given and let {6^)^ be a strict delta net with L^-bound K > 0. Define 
h := {x G M" : \x - Xo\ < b}, h := {x e W : \x - Xo\ < c + K\\F2\\i^,oo} and 
:= Ii X I2. Moreover set 



a := min 1, - — ■ — — EnrBli ' TiTrn • (^-2) 

/i,oo II -f^l II /3,0O / 

Then the initial value problem 

X'e = Fi{x^,X^) + F2{x^)6e, , . . 

Xe[-e) = xo, x^{-e) = Xq, 

has a unique solution x^ on := [— e, a — e] with Je), JJ) C J3. In particular, 
both and x^ are bounded, uniformly in e. 

Proof: We consider the closed subset X^ := {x^ G C°°{J^,W^) : x^{J^) C Ii,x^{J^) C 
I2} of the Banach space C^{J^,W^) with norm = ||a;||j5,oo + ||i||j5,oo- We define the 

operator : — )■ X^ by {t E Je) 

/t ps rt PS 

j Fi{x,{r),x,{r))drds+ J J F2ix,{r))5,{r)drds. 

First wc show that the operator A^ maps X^ to itself. Let x^ G X^ and t G Jg, then we 
have for the zero-order derivative of A^{xe) 

\A^{xe){t) - Xo\ 

/i p s p s 

j \F^{x,{r),x,{r))\drds + j j \F2{x,{r))\\5,{r)\drds 

<a\xQ\+a^ ll^l||/3,oo + a 11-^2 ||/i,oo||5e||Li 

< adiol + ||Fi||/3_oo + i^||i^2||/i,oo) < b, 
and for the first-order derivative 

rt 



Id. 



{A,{xe)){t)-x^\< j \F,(x,(r),x,(r))\dr + j \F2(x,{r))\\5,(r)\dr 

< tt ||-^l||/3,oo + ||-^2||7i,oo||(^6||l1 < C + K\\F2\\l^,oo- 
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At this point we claim that we can find a sequence of positive real numbers (a„)„>2 such 
that Yl'^=2 <3Sn < s-nd 

WA'^iXe) - A^(?/e)||ci(Je) < an \\x, - 2/e||ci(J.)- 

So let N 9 n > 2, Xe,ye ^ X^, and t G Je- Denoting by [n /*J the n-times iterated 
integral we obtain (with Lip(Fj,/j) a Lipschitz constant for Fj on Ij) 

\A^{x,)-A^{y,)\ 

<[2n [ ]|Fi(xe(r),x,(r))-Fi(y,(r),^,(r))|d2"r 



— e 

t 



+ [2n J ]\F2{x,{r)) - F2{y,{r)m{r)\d^^r 

< (Lip(Fi,/3) [2n y"V'V 
+ Lip(F2, h) ||5,|Ui [(2n-l) \\x, - 2/,||ci(j.) 

< (^Lip(Fi,/3)^ + Lip(F2,/i)||5,|Ui Jj^—^^)h^-y^\\c'{K 
Furthermore for the derivative of A" we find that 



<[(2n-l) j ]|Fi(a;,(r),i:,(r))-Fi(y,(r),y,(r))|d^"-^r 

+ [(2n-l) f m{x,{T)) - F2{y,{r))\\5,{r)\d^^-'T 
< (Lip{F^,h) [(2n-l) ]d' 



+ Up{F2,h)\\d,\\Li [(2n-2) J ]d2"-V)||x,-|/,||ci(j,) 

- (Lip(F„/3)^^^^ + Lip(F2,/0||5e|Ui (^^^jll^e-^ellcK.e)- 
Summing up we obtain 

||A;'(xe)-vi:(?/e)||ci(J.) <4max(^Lip(Fi,/3),i^Lip(F2,/i)j - 2/,||ci(j,), 

which proves our claim. 

Now Weissinger's fixed point theorem provides us with the existence of a unique 
solution Xe ^ X^. 

Finally, since x^. x^ stay in Ii respectively I2 (which arc defined independently of 
e), and x^ are bounded by b and c + i^||F2||/^^oo; respectively. □ 
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